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I. INTRODUCTION 


Hydrographic surveying includes many branches of science 
foc the purpose of the production of nautica T е ате 
specially designed for use by the mariner. The e r N 
tion of position in hydrographic surveying is as Importa 
as the measurement of depth. Before determiring an accurate 
Lydrographic position, accurate geodetic positrons sm 
shore control must be established. There are many methods 
available to establish geodetic control in the survey area. 
These methods are: 

Те Tr»xanguls tron 
2... r Trllateration 
3. Traverse 

4. Intersection 
5. Resection 

The process of making a proper nauticai chart consSisSmeg 
first of all, in setting up a Eramework of marks on the 
ground. Before 1950 the main framework of a geodetic survey 
almost always consisted of triangulation, which was repiaeed 
by traverse if the topography made triangulation impracti- 
cable. During the last decade, the introduction of elec- 
tronic distance measuring (EDM) equipaent has made both 
trilateration and traverse economical, and an acceptable 
substitute for triangulation. In fact, it appears prob min 
that these new methods will replace triangulation as the 
Main framework for new geodetic surveys [Ref. 1]. It is 
evident that triangulation and traverse are the main methods 


used for establishing control. There seems to be no 


1A position of a point on the surface of the Earth. 
Ec s in terms of Joos Ete eee geodetic longitude 
and geodetic height. A geodetic position impiles an adem 
geodetic datum. 
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Peeememt among the various agencics as to which of these 
two methods is mostly used. The U.S. National Ocear Service 
Pee ec Oo рс majority Of rts horizontal control surveys 
BOE ryarograpny with traverse (about 905) [Ref. 2]. The 
main factor for the selection of one or the other method 
depends on the geographical configuration of the survey 
project area and the availability of good EDM equipment. 
Traversing is a convenient, rapid method for establishing 
MeEazontal control. It is particularly useful in densely 
built areas, when» the coastline tends to be even, along a 
railroad track, and in heavily forested areas where lengths 
ENSTght are short so that neither triangulation nor 
trilateration is suitable. 

The objectives of this thesis are to show methods of 
Fa verse station position computation and a comparison of 
methods of least squares adjustment. Closed traverse 
Station positions were computed and adjusted using the 
Approximate Method and Ey the Least Squares üethod in the 
Universal Transverse Mercator (UTM) grid coordinates. The 
computer programs were written to calculate the adjusted 
traverse station positions. Test data included those data 
obtained during the Geodetic Survey Field Experience course 
Ere Naval Postgraduate School» (NPS)eemn October 1983. 

Bre sbesults OL Comparative computations are shown to 
more significant figures, in this thesis, than are normally 
considered desirable in production work. The same observed 
data are used with each of three computational methods. It 
is important to recognize that what is being compared is not 
observational precision but computational precision. Hence, 
it is considered necessary for a rigorous couparison of 
Semputational precision, including round off error, to show 
BeoeurtS to several more decimal places than is justifiable 


based on observational precision alone. 


TI 


A. GENERAL 


Traverse is a method of surveying in which a series of 
straight lines connect successive established points along 
the route of a survey. An angular measurement is taken 
using a theodolite at each point where the traverse changes 
direction. Distances along the iine between successive 
traverse points are determined by EDM equipment. The points 
defining the ends of the traverse lines are called turning 
points, traverse points, or traverse stations. Each 
straight section of a traverse is called a leg or a traverse 
line. 

A closed traverse originates at a point of known posi- 
tion and lis closed or arother point ol ktown HobNZzomu 
position. Traverse I-1—-2-3-F originates at pornt MW н 
backKsight along line IA of known azimuth and closes ono 
Fy. With a foresight along line FB also of known azimuth 


(Figure 2.1). This type of traverse is preferable to ali 
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Figure] Closed Traverse. 


T2 


ПЕЕ computational Checks are possible which allow 
Шс ошо sySteMatic errors? in Doth distance and 
direction. 

A closed-loop {closed polygon) traverse is a special 
case of a closed traverse in which the originating and 
terminating points are the same point with a known hori- 
Memeal position. Traverse I-1-2-3-4-I, originates ard 
terminates on point I (Figure 2.2). This type of traverse 
permits an internal check on the angles, but there is no 


check on the linear measurenents. Therefore, there is a 


e@estbility that an error proportional to distance may occu 


and not Ге detected. 


| 
| 





as АН АГАШ 


Figure 2.2 Closed-loop Traverse. 


Bee Oo NGLE AND DIRECTION MEASUREMENT 


Angles and directions may be defined by means of Ebear-. 


angS, azimuths, deflection angles, angles to the right, or 


interior angles. These quantities are said to be observed 


eSystematic errors may be caused NT instrumerts 
ог xctors SUCh as temperature or humidity changes which 
affect the performance of measuring instruments. 


ie 


— Const manca. UU... U U Ul... U... U U... u. U а uJ 


T^ 
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| 
| 
, 
| 
| 
| 
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when obtained directiy in the field and calculated when 
obtained indirectly by сори саак“ 

A theodolite 1S an instrument desijnec to observe hori- 
zontal directions and measure vertical angles. It consists 
of a telescope mounted to rotate vertically on a horizontal 
aXis supportel by a pair of vertical standards attached toma 
revolvabie circular plate containing a graduated circle for 
observing horizontal directions. Another graduated arcae 
attached to one standard so that vertical angles can be 
measured. 

Repeating and direction theodolites have features that 
are common to Doth types of instruments. Repeating theodol- 
ites are read directly to 20" or O1' and by estimation ME 
one-tenth the corresponding direct reading. Direction Tiik o 
dolites are usually read directly to 01" and can be esti- 
matec to tenths of seconds [Ref. 3, p. 215]. i genere 
direction theodolites are more precise than are repeating 
theodolites. 

The direction theodolite observes directions only and 
angles are computed by subtracting one direction from 
another. Assume that the horizontal angle AIB (Figure 22 
is to be measured with a direction theodolite. Тһе theo- 
dolite is set over point I, leveled and centered, anda 
sight is taken on point A. “he horizontal circ lens p 
viewed through the optical-viewing system and the circle 
reading is observed and recorded. Assume the reading is 
450 02' 40". The telescope is then sighted on point В. ШЕ 
horizontal circle is then viewed through the optical-viewing 
System and the circle reading is observed and records 
Assume the reading is 1249 11' 59". These two observations 
constitute directions which have a common reference direc- 
tion that is completely arbitrary. The clockwise порох 
angle is 

i = (124° 11° 59") - (459 02' ЧО = OoOO E 
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Ewgume 2.3 Horizontal Angle. 


Mee DISTANCE NEASUREMENTS 


There are several methods of determining distance, the 
Choice of which depends on the accuracy required anc the 
cost. For example, tacheometry, tapiny, and EDM equipment 
can be used. The general availability of EDM eguipment has 
practically eliminated the use of taping for the measurement 
of traverse lengths. Accuracies are comparable, or 
Superior, to those obtained with invar tapes. 

Distances measured using EDM ezguipment are subject to 
EUEOES arising from the instrumental components, calibration 
w te CCulpment, inaccuracies in the meteorological data, 
elevation discrepancies, and centering of the instruments or 
reflectors. The reduction of measured distances invoives 
converting the slope distance to a horizontal distance, 
Seqverting the chord distance to an equivalent arc distance, 
and reducing the arc distance to the ellipsoid. The reduc- 
tion of slope distance to horizontal distance is necessarv 
to compensate for the difference in elevation of the end 


points of the measured line. The horizontal distance is 
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reduced to the geodetic distarce Ly applyrnomu M 
corrector ard a chord-arc coEPErector to спел es 
distance [Ret. 8, pp qu 


D. ACCURACY 


In general, the accuracy of a traverse is judged on the 
basis of the resultant error of closure of the traverses 
This resultant closing error is a function o£ the accur G 
in the measurement of directions and distances. The classi- 
fication, standards o£ accuracy, and general soScurrc шы 
for horizontal control have been prepared by the Federal 
Geodetic Control Committee (FGCC) ;Ref. 5] and have been 
reviewed by the American Society of Civil Engineers, The 
American Congress on Surveying and Mapping, and the American 
Geophysical Union (Table I). “he third-order traverse class 
I and class II are of particular interest because these are 
the orders of accuracy the hydrographer is usually required 


to accomplish: 


Es. DATA ACOUISITION 


The data for this thesis were acquired at Moss Landing: 
California, by NPS Hydrographic students during the Geodetic 
Survey Field Experience course in October 1983.9 All Gaga 
known station positions and azimuths (Table II) were 
adjusted by the Coast and Geodetic Survey by third-order 
methods [Ref. 6]. 


The distances were observed by using a Kern DM102 and a 
Tellurometer MRA5. The Kern DM102 is an electro-optical 
distance meter. The measuring accuracy was +(5 om + 5 ppm). 
The Tellurometer MRA5 is a microwave instrument. Precision 
in terms of probable error, in the temperature range of 


-320C to +4¥u0C, is +{5 cm + 10075 pen) Stor aes EE 
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СОО Л ПЛ OF KION Stations ard AZinuths 


lc ae = lc pe eee OTS SS s AM 
| ШТ м Coordinates 
| Seat 1on ae m Bg dui Easting | 
| m. 
| 

| Moss 2 155, иви 608,279.04404 | 
| Ho In HO 2253531751 r o 00 005 65256 | 
| Grid azimuth clockwise | 
| From North mo gm tu | 
| P, | 
| From Moss 2 to Pipher stations BOUE 2327/78 | 
| Br Om Foln to Moran stations feces oe 262554 | 
P | 

——— YUY I PYP нта та mm hcic e] 


determination. The distances were observed in the field, 
corrected Бу temperature and pressure for propagation error. 
Unfortunately, meteorological data are generally accuired 
only at the end points of a measured line. Using the mean 
of these meteorological values only approximates the actual 
conditions of the entire measured line ani does not 
Enpletely correct for errors in the propagation velocity of 
electromagnetic radiation. By applying the elevation, sea 
ievel, ard scale factor corrections they were reduced to UTM 
bea distances [Ref. ü, pp. 124-125]. The UTM grid 
distances and standard deviations of the distances were 


determined in meters (Table III). 


The angles were observed by using a Wild T-2 theodolite. 
Mesure the Correctness of the begiarirgz and ending 
azimuths, a check azimuth to a second station of known posi- 
tion was observed. The angles were observed at stations 
Pees 2, HMHossback, Dune Temp, and Holm (Table IV). All 
observations were made by NPS students and conform to 


Cci fication for a third-order class I traverse. 


iS 


TABLE ETI 


Grid Distances and Standard Deviations 


{0 q TE ae aue aM OM = n — :  -— ee ee ee ee ee ЭЖЭ 


Grid distances Derea 


Mossback 


Dune Temp 


Stations 


and 
and 


and 


Mossback 
Dune Temp 
Holn 


m 

| 

| Woss 2 
| [ 

| 

= 


The 


г ——— 


| 
ПЕЗЕКЕТа ВЕ 
station 





Pipher 
Moss 2 
|Mossback 


Dune Temp 


' 
Moss 2 

Mossback 
Dune 
Holm 


Tene 


| 

| 

| Standard 
Distances deviatior 

( m. ШЕ | 

1,2 0. 00m | 

30527743 0.00 | 

6,0 769 m 0.008 | 

— — — — € M — 


TABLE IV 


Observed Angles and Standard Deviations 


a Se me ee rm m ar 


Foresight 
station 
Mossback 
Dune 
Holm 


Moran 


Temp 
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Observed 
angles 
О t " 
246 05 3252000 
2з 08.600 
190. 15 023600 
270 NEEDS 17291900 


Stand ш е 
deviation 


3 
| 
! 
| 
| 
" 
01. 984 
01.495 


206 


| 
| 
| 
| 
01.6185 | 


Meee SE CONPU TATIONS AND ADJUSTMENTS 
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À. Pera TRAVERSE, COMPUTATIONS 


О Ел уе кс ОЕТ ае i) Obiginates at Statior 1 (known 
position) and terminates at station 4 (also known position) 
ELS. —5 compute the forward azimutk (Table VI) of an 
unknown leg, the angle is added to the back azinuth of tne 


previous leg (Zzquation 3.1). 


! 
Az, = Az. + c EE (ne 
Where: i = the number of legs, 
Az, = the forward azinuth of an unkrown leg, and 
Az. = Mien ео ааста 1 azimuth. 


For the computation of UTM grid coordinates, let QE; 
and Q\N, be designated as the departure and latitude for leg 
i (Figure 3.1) so that the general formulas to compute the 


departure and latitude are 


QE, E EIE (3:2) 
AN, = di COS А2. (2573) 
where x rox MSIE EE C ME. 


n the number of the observed angles. 

The algebraic signs of the departure and latitude for a 
traverse leg depend on the signs of the sine and cosine of 
the azimuth of that leg. The algebraic sign of the depar- 
ture and iatitude is determined by the following rules: 

PSR UIN grid azimuths are refered to grid north. 

2. For azimuths between 0° and 1809, the departure is 
pius; for all other azimuths, the departure is minus. 

3. FOr azimuths between 90° and 2709, the latitude is 


ДЕШЕ ОЗЕ оет ага оше, Ее latitude is plus. 


zl 
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| 

ei 
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Figure 3.1 A Closed Traverse. 


The calculators and computers with internal routines for 
trigonometric functions yield the proper sign autonar Te TES 
given the azimuth. The major portion of traverse conp Ti 
tions consists of calculating departures and laticu M 
successive legs and cumulating the values to determine the 


coordinates for consecutive traversen NN DEM 
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TABLE V 


Data for Initial Traverse Computations 


SS <= Te? eS em -—— — s — . r F  — — — —- 


UTM grid Northing J, oU 5505550 5 7 06 m. 


E 7 | 
| 
| Leg Distance Angle | 
m. 
| O 1 " | 
1-2 d, = 1,424.004 OU = 286 65 83. 200 | 
| PIC ce = 365.744 aA, = 22? m03. 600 | 
| 3 = 8 da = 6,%76.271 ©- = 190 15 02.600 | 
| 
| 
| ae initial UTM grid azimuth Aze = CORR T6 23. 178 | 
| 
Station 1 | | | 
Шиа са каѕсіда E, = 608,279. 01449054 в. | 
P EE  — 
The coordinates of station j (j = ií + 1) are E; ana 
Ni the coordinates of station j are 
DLE CC (3.4) 
Б Е: ÔN, (35) 
Mieke stea 1 72,3, + - » i. 


The values of departures, latitudes, and coordinates are 
Shown in Table VII and were computed from the data in Tables 


ШОШО Т OY application of Equations (3.2) through (3.5). 


Eee COMPUTATION OF DISCREPANCIES 


ОШ ose raverse at Noss Landing (Figure 3.1) origi- 
F Led at station 1 (Moss 2) ard closed at station 4 (Hola). 
The closing errors for a traverse are caused by observation 
errors in the observed angles and the measured distances. 
Ir Closing errors may be computed by applying Squations 
ER through (3.6). 


2:3 


TABLE VI 


The AZinuth Са сита ы 


[o UU те —————————— UTERE 9 Ъъһ" 

| rorward Back Observed | 

| Azimuth Azimuth Angle 5 | N | 
Az. 100 16 25 EI: 

| б, 246 05 33.2900 | 

| Де, 346 22 06.978 | 

| - 180 00 00.000 | 

КЮ 166 22 06.978 | 

©, 222 51 08.600 | 

А2, 389 13 15.578 | 

- 180 00 909919 | 

| Iz 209 13 15.578 | 

| a, 190 15 02.600 | 

| А2. 399 28 18.178 | 

| - 180 00 00.000 | 

| лса 219 28 есе | 

a 


The angular error of closure may be computed by applying 


Eguation (3.1) becomes 
n 
Az. * 2 O- (n-1)1809 - Az, = ñ, (3.6) 
i -1 


where AZ, 15 the fixed closing aZimuth апо и ие 


angular error of closure. AZ, is the grid aziuuth from 


stations Holm to Moran (AZ, = 1369 33' 26.334") (тареш 


The computed closing azimuth 1s equal to А2 43 plus the 
observed angle at station Holm = (2199 28' 18.178") + 

(2779 O5' 17" = 4969 33! 35.178" = 136° 33° 3350 7 P АН 
Equation (3.6), the angular еггог о ci ucl 

(1369 33' 35.178") - (136° 33% 26.5338") = +3023 ee 
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ОЛООР 


Initial Traverse Computations 


EEUU a -————— 
| Values in meters Values in meters 
| N, uomo oo) 20.6 E, 608,279.04404 
| SN. 13e 92 OR ж J 60156 
| N, 1507 a773 E, 607,943. 44238 
AN, J 5700071 AE, 5515571 
| N. 4,074, 258.94534 ' ЕЗ 608, 121.99 109 
AN 4,999.28284 А Ба ПИБ - 94 255 
| N. WO 2539.22818 Е 6127239. 955515 
| a Sree 
| 
a L «EE ERES m —— 
Coe on ute the errors in closure in position for a 
eosed traverse. Equations (3.2) to (3.5) are applied to 
brain: 
n -1 
me en’ (3. 7) 
n-1 
or | 21 
Wo = E, Е (3. 7а) 
н + SON ME N- (са) 
Nhere: n - the number of observed angles, 
ома Decal culated discrepancies, 
PSU - е хез Еа 1 собЕа тіпа гес; 
ЕТ, N. = the fixed closing coordinates, and 
БОШ ОН = cre CONDuUtCd=closing coordinates, 


E result of the errors in closure in position of a 


meaverse is 


25 


ba a ee eS a >>> wa — ай 


Cs d = KS + HS ) 142 (3. 9) 


The computed closing coordinates of this traverseucme 
Bac 612,238.93864 m and N, = 4,079,258.22818 a. The fixed 
closing coordinates are E, - 612,238.85256 m and N. - 
4,079,258.31754 m. The precalculate4d discrepar ii- cE 
computed by using Equations (3.7a) and (3.8a). The results 
0.08608 m and W, - -0.08236 m. By applying 
Equation (3.9), /MAd is 0.124008 m. From Table Iii, the tome 
measured distance is 8,266.019 n. The ratio of the distance 


error of closure to the total distance is 0.124083/78, 26670 WT 


Il 


are Wo 


or 1 part in 66,618. The ratio is ад шл сак о ы eee 


goodness of this traverse. 


C. ADJUSTMENT OF A TRAVERSE BY AN APPROXIMATE Sates 


The angular error of closure of traverse may be distrib- 
ute! equally among the observed anyles. The angular error 
OF closure of this traverse is +8.344", which corresvonds to 
a correction for each angle of -2.211'". The observed anaes 
were correctel by this value. The adjusted azimuth was then 
computed. The departures, latitudes, and coordinates were 
recomputed Fy using the adjusted azimuth (Table VIII) MRE 
Closure corrections in E and N coordinates are €0.096559 
and -0.04326 m, respectively. The resultant closure ще 
0.10562 m. The ratio of the distance error of closur MEN 
the total distance is 0.10562/8,266.019 or 1| part In SPP EMEN 

The adjustments of a closed traverse by the Approximate 
Method are completed by the compass rule which proportions 
the errors in E and N coordinates according to the dus 
of the course [Ref. 4, p. 354]. The corrections are yao pies 
to the departures and latitudes prior to Como tr 


Саса ес 
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E d 


ТАЕ ОЕ 


Initial Traverse Computations. 


AN, 


= 
| 
| 
| 
| 
| 
| 
| 
| 
| 


AN, 


CN. 


| 
| 
d 
| 
d 
| 





Adjusted 12M 


Observed angle Correction 
AZo 
ОШОО 52" 2.21] 
AZ, 
A Z5. 
ООО КОТ ӨЫ" — 2,211" 
AZo 
А2-_, 
ОООО ИОВ и — 2,7 ү" 
AZ, 
anes 
Poe Ces 17.0" = 2,211" 
n 
Values in meters 
N, ОЛ 2 55785206 Se 
Ss 07 ae 
N, i dg MOON TEES 2, 
2V9 2044 E 
N. ШО ЛШ Б O1 6657 ЁЗ 
4,999.41523 aver 
Ny ad 59. 35080 E 
N- 0791799231754 Er 
dn C203 26 dE 


5 Method with 


Adjusted azimuth 


O q t 
MODE I6» 3.778 
206 05 40.989 
346 22 94.767 
- 180 00 00.000 
ТООМ o OR 767 
OPNS OO. Ou 
EOS MEC 11.15565 
*—190 00 00.000 
DOSES 11-159 
19:0 015759091083 89 
3950 29 —11.5455 
- 180 00 00.000 
DUCI DESI SI5 
2T EDS EDI 789 
o Өе CO SH 
Values in meters 
9B 279 901): 018 
3355261659 
SOT SU Jun 
178.9594187 


6005121 06512 

(1716578675 
Gale 3 Oar) 56 21 
GNE 992590 


a 
pO Ce. 


an  — — — ——— 
= oS —— rss N — s ee SS aa cs ee ee EIN EEE EES =. 


bo mara a amra amra аа шз ше es, o E ee ee es, ыр eee Sonate es ED N... s... --. ,. p, TOG r, OGD „ч s r OG D, a... D Ds LL a D cee SE cere meee — — _ — аә ES 


Ô E, = ( аЕ 0587 (3.10) 

ON; = ( dN /D) . d, (3.11) 
Where: OE, = correction to a x 

ON, = соггес осо SN: , 

dE = total closure correction in the = coordi mame 

dN = total closure correction in the N coordinate, 

and 
D = total distance. 


The adjusted E and N coordinates (Table IX) were 
different frou the fixed E and N coordinates (Table un D 
0.00001 m due to round off error. The calculations 
adjustments were illustrated in sections À, B, and C by 
using the hand calculator (TI-59) and rounded oom uu 
decimal places. The computer program was written to calcu- 
late the adjusted traverse station positions by the 
Approximate Method by using 16 decimal places (Appendix A). 
The approximate traverse adjustment 1s based on an assumed 
condition that the angular precision equals the precision in 


linear distance. 


D. ADJUSTMENT OF A TRAVERSE BY LEAST SQUARES METHODS 


The method of least squares provides a rigorous 
adjustment and best estimates for positions of all traverse 
Stations. The Least Squares Method is used to 
Simultaneously eliminate closing errors in azimuths and 


coordinates of traverses. 
1. Ihe Principle of елы е 


The fundamental condition of the least squares tech- 


nique in surveying requires that the sum of the squares of 
the residuals be minimized. A residual is defined as the 


difference between the true and observed values. In making 
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TABLE IX 
Adjusted Coordinates (Approximate Method) 


m — eRe ——— B шшш зш e — ш шша сш ш є_ 


—— n — 


= 
| | 
| Departure Correction Grid Values in meters | 
| E, 608,279.04404 | 
| (SE,  -335.61689 * 0.01659 5535259590 | 
| E 607,943. 44414 | 
| QE. mo sus 0.000295 178.54613 | 
| Е. 600 212990 27 | 
DE, o e ЫЛ *00:97549 6116756228 | 

| Е 4 10. 238.09 285)9 | 
| | 
| Latitude Correction Grid Values in meters | 
| N, TO SS reao | 
| 8 11153832. 998907 = 10. 00745 ПЕ ОСЕ ЕС | 
| N 07505 NE 668 | 
| AN, 319.20844 - 0.00191 j 9270252 | 

N. n 0212595959 362] 

| AN R2 = 003339 15 9992558138 | 
| N. OTI 2555317255 | 
L | 
ЕЕ TT I LLL LL L. 22 


physical measurenents, the true values can never be deter- 
Mined. Тһе least squares principie establishes a criterion 
for obtaining the best estimates of the true values. 

If the best estimates of the true values are stated 
by x, and observed values by х., the residuals are expressed 
as 

v. M = x. 

The fundamental condition of least squares, for 
uncorrelated observations with equal precision are expressed 


as 
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n 
206,0? 7 (9)? * (9)? * (v,)2 + . . - + (P i A 
t= 1 
орпа ока отн 
viv = Minimum (3. 125 


In general, the observed values are of unegual 
precision. The observed value of high precision has a snall 
variance. Conversely, a low precision of the observed value 
has a large variance. Since the value of the variance goes 
in opposite directior to that of the precision, the observa- 
tion is assigned a value called weight corresponding to a 
quality that is inversely proportional to the observatior's 
variance. 

For uncorrelated measurements Хр, £54 X36 2 + © 
„œ; respectively mimm 
weights ог these uncorrelated measurements are, 


2 2 2 2 2 2 
P E oy © # EE Ow o Pas = Oy Ó., ° ° . • 


11 


| | 2 2 2 
X4, With variances Се б, б, scusa MS NO 


2 2 
DICA (3-98) 


2 
where O, Ís the proportionally constant of an observatis 
of unit weight [ Ref. 7, p. 67]. These weights may be 
collected ints a corresponding diagonal P matrix, called the 


weight matrix: 


or Ota Е 0 
Ti 
0 0 : 0 
0 y 0 
P437 7 
P = е е е ə е a 


nn 
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The weight matrix consists of weighted observations 
Bp traverse with mixed kinds of measurements. They are 
Gistances and angles. The variance of a measured distance 
is expressed in meters and an angle in radians. Гог 
meoir elated observations of unequal precision, the funda- 


mental condition of least squares is expressed as 


* 2 2 + 2 + ye * S ° d Y + v2 
E. Pu * Pour Р"; P33 3 Pm n 


Minimum 
meetin matrix forum 
vipPV - Minimum (3.14) 


In the problems involving the adjustment of observed 
values, all observed quantities are expressed by functiors 
Of the quantities to be determined. In simple cases these 
relatiors are linear, but when this is not the case, the 
relations must be converted into the linear form by 
expanding them into Taylor's series. The terms of higher 
‘order are neglected, so as to obtain linear relations, 
solving the resulting linear equations, then iteratirg until 
the effect of the neglected higher order terms are 


Meni mized. 


2. Least Squares Adjustment of Indirect Observations 


——— — сс A —— — — = с шс — —— m ë - шы шр ш == mmm mm mm f m Í чш == чш <= шшр сэт Q wm RED 


Ihe observed values are related to the desired 
unknown values through formulas or functions which are 
called observation equations. One observation equation is 
written for each measurement. To solve for the best value 
of each unknown parameter, at least one redundant observa- 
tion equation must be written. That is, the number observa- 
tions must be greater than the number of unknowns. The 
linear observation equations can be written in the general 


fern as follows: 


2 


TO зи V2 m р ТЕ 1 1 
m ° ° - + = + 

EU a55X5 f ат Хт 22 5 80 

a ot с^ ? алоо Т - i акы = Gu Е Y. 


where n is the number of observations; &@ is the number of 


the unknowns; ас; axo: + * © āno are constants; a4, 8415, 
- . . 4pm аге coefficients ofiti-mun NL LO X, E 
X4; and Vj, V5, + + + Vy are the residuals. 
Because the observations G, (1 = 1, 2, ОПЕ 


not free from random errors, each G. must be corrected by a 


residual value, у. let b. SECO. а; 9: Thus 
аут ху + арх +... + ај Хн b, 7 Vj 
аху + арх +... ваха В = У 
+ + + = = 
а пт 1 ANDES 57—65: ES m b. M 
Or in the matrix fora 
V - AX- B (3.15) 


This equation is called the observation equation or 
the observation equation matrix. 

For uncorrelated observations o£ unequal precision, 
substitute the value for the V matrix from the observation 


Equation (3.15) “into Equa Vong NN 


В 


PV - (AX - B) P(AX - B) 


[ (Ах)! - В ]Р(АХ - B) 


(WI AT —- Bl)P(AX - B) 


(ITM p – ВТР) (АХ - B) 


mbi »pe- plpAX - BPB 


EMEN i pp = xr A PB + B'PB 


X M PAX - 2x'alpp + Bl pp (3. 16) 


(from matrix algebra, (AX)! = x'a! anā BĪPAX = X'A' PB] 

The minimum of this function can be found by taking 
the partial derivatives of the function with respect to each 
unknown variable (i.e., with respect to the X. , Xər * + > 
eee must equal zero. Hence, 

O (vipv) = 2АТрАХ - 2a7PB = 0 (3.17) 
Ох 

Puvuduneeseusthon (3.17) by 2, the following result 

ms Obtained: 


T 


A PAX = A’ PB 


This represents the normal equations, and by multi- 


plying this equation by (aA pA) -1, the solution is obtained: 
X = (a!'pay—1a'pp (3.18) 


For uncorrelated observations with equal precision, 
meen welgnt matrix i5 an identity matrix, and eguation (3.18) 


becomes 


X = (АГА) -1АТВ (3. 19) 


This equation can be derived similarlv to the 


K (al precision case. Equations (3.18) and (3.19) are the 


Lasic least squares matrix equations. 
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If the relations are nonlinear, the сеа опы иш 
be converted into the linear form by using Taylor's series 


expansion [Ref. 3,9p^ 99079] MEM (X,, Хх, +. x 


be the general observation equation that is a nonlinear 
function. The Taylor's series expansion is 


pes £(x,, ©) 7 n ғат ох осор . o. + Of Ox 
ox O Qx? 





+ Higher order terms 


° ° 


ki re к; X», + « « X, are the approximate values oc usm 
variables at which the function is evaluated as os: 
traverse problems, an approximate value can be the 
precalculated value, and Ox,, One СЯ: DX can be the 
corrections. The higher order terms in the series are 
neglected and only the zero and first order terms are 
maintained, the approximate value must be improved by 
successive iterations until the erfect of the neglected 
higher order terms is minimized. After linearization eee 


observation equations become: 


v = Lr X5, gae x° ) +O£,5x,+ O£,dx,¢ S + Of o G, 
0x 0x, O 

y Eo. X», 0. XA) * 9£)0x, « IE dxat а: 2220.7 G, 
Ох, OP Х т 

yen RC X^, . хач) * S£, 0x, « BE dxat a ЖЕ * of bxc Go 
Ох, Ox", x= 

OL ln the natrix Form 

V = AX-B | (3. 20) 
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where 























fone tu ТРОЕ V 
1 l 

2 х% х, Ox? 

Of DEMENS. . . . 9f y 
2 2 

0 x* до, Ox 

A = Е 2 а V = ° 
Ое сое а: © É. Vn 
Ox, д X ox 


G, ы Í, (x , X.” = . ә ° Kan ) $ i. 

б = Ё. (х, Xor ee. rk.) Ox, 
B = е E е 3 ° X = - 

CB Mex Lor eee 4 Or 


The remainder of the least squares procedure is the 
Беше а5 indicated by Equations (3.16) to (3.19). 

There are two types ots observation equations for the 
adjustment of the traverse.  Thev are the angle and distance 
observation eguations. Both of them are nonlinear, and they 
must be linearized using the Taylor's series. 

To derive the observation equation for the angle ©, 
EN ure 3.1), first write the two azimath observation 


equations for the azinuth À Zik and AZ; ,Q where AZ.) is the 


Buck azimuth and Az, а= спе ягонага azimuth. 
b E, — Ej 
AZ... — arc л 
E: INE 
AZ.. = 


| arc бап р 
) i 
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Then, the angle observation equation for пешашон и 


V. - AT Saat ACC > e. 
| l] I K I 
Ej- Ei E NM а 
- arc tan = - arc tan = = (3-290 


Equation (3.21) is nonlinear in the parameters. 


Write “this equation Gn guns. iounrorm 


Thus, the linearized form of the angle observation 


equation is 








v, = f. (E, N, E., Nr, BE, N. ) + Gf. 05. СОО 
I l I I J ) K k zi I = | 
дЕ" дм 
+ @ е. бЕ. + Of, ON, cba lem + Of, ON, ECCO 
дЕ; | ONS дЕ" ON’, 


The observation equation for the distance d; between 


two points. i amd] = 


ү, = ( (E. - ESL (N; L ` 


OL ir the functionali sem 


Thus, the linearized form of the distance 


observation equation is 


< 
II 
rh 
tj 
o 
- 
д 
о 
= 
td 
o 
Z 
e o 
ee 
+ 
Q) 
Fh 
C») 
[rj 
+ 


ОТО 
ON. 





ES 
с) 
m 

=e О 


OEE O = a. 
J Эм J l 
j 





For the adjustment of the traverse which was 
conducted at Moss Landing with the technigue of adjustment 


of indirect observations, four angle and three distance 
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observation equations need to be written. The sever equa- 
Ioue vouenelude ас unknown parameters the four coorii- 
Mates Oc Stations 2 (Mossback) and 3 (Dune Temp). Four 
Шона ecuations are formed and solved for corrections to 
the approximate values (precalculated vaiues) of the parane- 
ters. The corrections are added to tne approximations to 
update their values, and the solution is repeated until the 
last set of corrections is insignificantly small. 


the linearized eguations of this problems are 


т, = з + a,, ON, fiae OE, + а,, бн. - b, 
= Tt = 

V5 a, OE, + a, ON, + а,. D + a, , ON, b. 

v, = ОЕ. + a. ON, + as. Sis + à, DN. - b- 


eeo der cerns oE tho angle and distarce 


functions are 


Е-Е 
= к= Ше 
EE arc tan =n Az- 
pe re u o 
f. = arc tan B c = a tease an E 2 
2 N$ - N3 N1— N5 
E B Ес, E' 
f, = arc tan N = -are am 2—2 
Í ES— E i 
f - oa oe 
f, AZ arc tan TE 
ЕЕ Ер (N = в ) 24172 
LOS Е 02 + ELI үз: 
LIU LO Ee Емо 21722 


By using the data of the known coordinates (Table II) and 


the approximate values (Table VII) 


а, 9f, - Nó - N = 0.68246 x 10-3 
дЕ? ° 2 ° 2 
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ло 
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= or ss E2 - E; = 0.16550 ЕКШЕ 
END (м5-№)2 + (Е5— Е, ) 
= © Ё, — 
OE 
= Of, = Q 
ON" 
=0#, = - N3 - Nb + N, - М5 | 
Е (N3- N5)* + (Е5- Е, ° (м, №) (Е, є®)° 


2 3.00569 0019. 


deo L- 422.2 p E 
ә ° ° o о 2 E ° 2 _ =° 2 
2 [wN$- N20? + (E2- E.) му №,)2 + (Е, Е°,) 


= 1.16926 х 10-3 


2.38062 1990 


l! 





= % f = | N^ = N^ j 
( 














Cee сме ес зеш д 
ЁЗ к = 
=?f, = f E3 - Eo. i = 1.33876 x Wee 
° ° о \2 o о 
ON? (s EE (CENE 
RT Nd — N | 2:38621 x 10-3 
DEN (N5 -N3 + (E? — £^? 
=Q £f, = - E2 - Е“, |` - 1.33476 0] 2€ 
° ме 32 ° =° 52 
BN (м5-№ )2 + (Е5- Е5) 
al N4 - N3 l| N5 - Ny | 
e 2 N$ 32 I EU S cM EE 9. ONE 
3 (Ма №3) + (Е, ЕЗ) (N. No + Eo ee 
= =) 2290545 1 sane 
k | E4 - ЕЗ l- E5 - E^ | 
o N° — a EU Е ТА MEC о ° 2 
3 ЕМ CN ME смо, )2 ЕЕ а 


= 1.43291 х 10-3 
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Q? 
rh 





Q 
zZ 


Q) 
Fh 





CD 
m 


OQ) ° 


il 


Q) 
lh 


°| 
ZA 
U 95 


© 


о 


= f N3 - N4 | 0.11920 x 10-3 
PME 2 QM e 
(eee ЕЕ, 


r E - Ey | - 0.09816 x 10-3 
2 
(N° ИЕ *(ES- E, ) 


EE 550066 x 10 = 


EY = ED : 
(E,- E HW ЕТ БШ г 


тт - N2 = 971.83201 x 10-3 
ke, - E2)* HON, E 
0 
0 
| E> - ‚т - 488.17948 x 10-3 
o ° 2 Is 
(E5-E шм, 
№ - EHE 2с х 10-3 
— Piu )2 11/2 
- S ES |: 488.17948 x 10-3 
(E5- E39 e 5 - N32]! 7? 


“| No - N3 | 872.7%326 х 10-з 
° 912 ° QU OUT 2 


0 
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= 227 = 7 
№ о 
3 (Es — 


OF in the natrix roci 


"4 


| = Е) 


E3 - E4 


N3 - N3 


2 9 2 7A 
ъем | 


° 2 Z 
“Q eee 


0.68246 0.169350 0.00000 
-3.06868 1.16926 2.360621] 
2236621 -1.33476 -2.50541 
A = 10-3 0.00000 0.00000 0.11920 
-235.67466 971.83201 0.00000 
-4 48. 17948 -872.74326 488.17948 
0.00000 0.00000 -635.68254 
The compütarionmmer PETC S 
Irom Tables II; IIIS IV; anday 
b, = A, = Е. = 0.00000 radians 
b, = a, - Ё. = 0.00000 radians 
b. = QA- f, = 0.00002 radians 
b, = A,- £, = 0.00002 radians 
b. = 3, - £4, = 0.00000 meters 
De =u fec = 0.00000 meters 
b, = d - £, = -0.01426 meters 
or ir the matrix forn 
в = го о 0.00002 0.00002 о о 


= 635.68250 x UNE 


- 771.95059 x 10-3 


0. 00000 
-1.33476 
1.43291 
-0.098 16 
0. 00000 
872.74326 
-771.95059 


by uSing the daca 


-0.01426 i 


The weight matrix can be computed by using the vari- 


ance of the observed angles and distances from Tables III 


and IV by applying roseus mer, 


(oa 


observations of unequal precision, 


defined as the diagonal P matrix. 
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For uncorrelated 


the weight matrix is 


p, - 1/ sin?(1.988") = 10,808,537,426.79957 
е^ Smee ИЕККЕ ИО 1 552 70987 219. 02074 
А sin2(1.203") = 29,398,082,997.43487 
Е ЕРИСИ ee =) 16,332, 144, 194 . 087 30 
D 1 / (0.001)? - 1,000,000.00000 
p О 09001) = B 1,000,000.00000 
pee ee (0003) 2 = 111,111. 11111 


For uncorrelated observations of unequal precision, 
the values of A, B, and P matrices are applied to Equation 


INS. The correction vector is found to be 


OE, = 0.01284 
x= OR, = 0.00336 
SE, - 0.01079 
ON, = 0.00495 


The approximation values were improved bv adding the 
Berrections to the first approximation values. The improved 


approximation values after the first iteration are 


Eo = 607,943.44238 + 0.01284 = 607,943.45522 nm. 
М, = 4,073,939.74473 + 0.00336 = 4,073,939.74809 a. 
Еу = 608,121.99110 + 0.01079 = 608, 122.00189 п. 
N, = 4,074,258.94534 + 0.00495 = 4,074,258.95029 n. 


Using these values, the solution is iterated. After the 
second iteration, the correction vector is zero to six 
decimal places, so the improved approximation values are 
the final estimates of the coordinates. 

For uncorrelated observations or ecual precision, 
fre, COrrection vector is found by solving Equation (3.19). 


BMENM o orrectron' vector after the first iteration is 


GE =) 0.01756 
X = ON, = 0.00426 
SE, = 0.01660 
ON, = 0.00479 


4 1 


The imprcved approximation values aiter the SiS CER e 


are 
E, = 607,943.44238 + 0.01756 = 607,943.45994 1. 
N. = 4,073,939.74473 + 0.00426 = 4,073,939.74899 n. 
Es = 608,121.99110 + 0.01660 = 608, 122.00770 m. 
N, = 4,074,258.94534 + 0.00479 = 4,074,258.95013 m. 


The solution is iterated Fy using these adjusted 
values. The correction vector 1s zero to six decital рсе 
after the second iteration. These values are the final 
estimates of the coordinates. 

The standard deviation cf an observation which has 
unit weight can be found by the following equations, for 


uncorrelated observations of unequal precision, 


С = [RS | 1/2 (3:290 
n-m 
ror uncorrelated observations of equal precision, 
O, = Las 1/2 (3.23) 
n -m 


where n is the number of observation equations and m is the 
number of unknowns (Ref. 7, p. 249]. After calculatinmcMM 
best estimate values of the unknowns, the V matrix can be 
computed from Equation (3.20). 

The standard deviations of the best estimate values 


for the unknowns are then given by the following equations: 


GC s; 1/2 (3. 208 


| 9 


where б. is the standard deviation of the ith adjusted 
quantity. The cuantity in the ith Drow of (the unc NE S: 
is an element of the (AT pA) -1 matrix Сог ипсогге аре 
observations of unequal precision. For uncorrelated 
observations of equal precision, S: is an element of the 


ll 
(АТА)-1 паёгіх С Кее. 7 р. 22-086 
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"Es oS 
321 322 СГ oos Som 
(АГРА) 1 ог (АГА) –1ї = : | е, ЖО ОИ. 


m1 


The standard deviations of adjusted angles and 


ЕЕ са псеѕ сал be computed by the following equation: 


6. - б,| 0) 1/2 (3.25) 


Dre © is the Standard deviations of the ith adjusted 


Emu E2tyv. The quantity ln the ith row of the V matrix, Bre 


is an element of А (АГ РА) -1А7 matrix foL urcorrelated 
observations of unequal precision. For uncorrelated 


observations of equal precision, Q. is an element of the 


Ma armia matrix [Ref. 3, p. 9121. 


2n 
A (A! PA) -1АТ x š mu... 
Or = е ° О о ео о е 
А (АТА) –1А7 : д ле: е. 


пі п2 е = е 3 е О 


À closed traverse was conducted at Moss Landing. 
eect the second iteration, the V matrix can be found by 
BEEp Ong the values of A, B, and X matrices to Zquation 
E20). 
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For uncorrelated observations of unegual precision, 


ths 7V natrix is оос ое 


Y S 0.93109 x 10-5 radians or 1.92134" 

Voy Wes -1.63120 х 10-5 radians or 325990051 

v4 7 -1.28363 х 10-5 radians or -2.64768" 
VvV = V ES -2.30436 x 10-5 radians or -5.753059 

Y 75 0.00024 meters 

VES 0.00039 meters 

ie 0.00358 meters 


The standard deviation of unit weight is found by solving 


Equation (3.22)) thewnestite ac 


б. = # 2:69685 


The values of (A! pA)-1 are 


1.295 178. 0.22019. 717225959 5 OL UNES 
10-5; 0.22019 0.013919 0 2151 i КЕРЕ 
1.25575 0.21810 1.86H8081 10010 3 
0.10788 0.11127 00.03315 f0 o Ds lD 


The standard deviations of the best estimate values for 
positions can be found by solving Equation (35201 aw. 


results are 


a = + 0.00954 meters 
Ó N2 = + 0.00317 meters 
без = ж соо л TN 
o = + 0.00405 meters 
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The values of A (A PA) -1AT are 


0.064 
me 0 2 
2029.21 

0.010 
=. 165 
ZOG 
559 


19-10 


zou 


О 25 


=0 5 


02007 
3.415 
4.690 
41.59 


021 
UE 
ОЕТ И 
=O. 0 14 
Оба 


Og TO 
Dm 
ES) 
E 2 
21.913 
Ends 301 
ОЛ 99501 


bs 15m 
Ee 
2E EE 
= 
9917 
zm 
-1047 


xo 65.58 
g 5090 741.59 
s 7 17/9 

ESI =9.501 

т -10857 
О = 1525 

E925 ose 


The standard jieviations of adjusted angles and distances are 


Bounc by solving Equation (3.25), 


Os, = + 
Oro = + 
б.з = + 
O.4 = ¢ 
Ó = + 
Cio = + 
баз = + 


1.40346 
2.64103 
2379952609 
0726137 
000269 
0.00267 
0.00745 


Seconds 


seconds 


seconds 


seconds 


neters 


neters 


meters 


the results are 


For uncorrelated observations of equal precision, 


mre) vy Matrix is found to be 


ха 


waa 


v 
d1 

v 
d2 

3 


v 


а. 


ile 
БЕ 
= 
= 
O. 
Or 
0. 


26896 
56893 
75414 
23358 
14037 
23844 
24365 х 


м м K K Ki 


10-5 
10-5 
10-5 
10-5 
тот 
10-7 
То 7 


radians 
radians 
radians 
radians 
meters 


meters 


meters ` 


OL 


DE 


OL 


Or 


O 1712" 
99 20 Ei 
See oly 
-4.60710" 


uc -tandard deviation of unit weight is found by 


Solving Equation 


(351235 


РО у 10—5 


the result is 
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The values of (al ay~2 are 


0.28193 0.06827 -0.854 4 Spo n n 

1035 0. 06827 0.01664 -0.20731000 r ШЕ 
-0.85449 -0.20731 2.59561 EE 
0.70380 0.17075 -2.13/25 EES 


Tha standard deviations of the best estimate values for 
positions can be found by solving Equation (3.24), the 


results are 


O-, = + 0.00107 meters 
Oyo = + 0.00026 meters 
бе = + 0.00324 meters 
Олз = + 0.00267 meters 


The values of А (АГ А) –1АТ аге 


1.472 -27.197 26.960 -1.2395 0.5 50-0 SELON 

-27.197 503.1 -498.7 22.855 3.3408. 0 150 NE ECL 

26.960 -498.7 494.4 -22.66 0.5700 05523 20 ESL 

Шз “1.235 22.855 -22.66 1-038 0229001 0S E E 
=0 51 0.348 0.370 -0.201 10000 -O 00 UM 

20.093 0.480 0.523 -0.159 -0.001 оси E 

2051856 0.573 9.518 -0.134 -0.001 DUDEN Е 


Ihe standard deviations of adjusted angles and distances are 


found by solving Equation (3.29), theres h oE E 


O4 = + 0.15917 seconds 
Oi = + 2.94270 seconds 
б.- = + 2.91732 seconds 
oF = + 0.13370 seconds 
Ou = + 0.00002 meters 
on = + 0.00002 meters 
or = + 0.00002 meters 
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The computations in this subsection were performed 
EE -SNPSOEBM 30339eo0mputer using 16 decimal places and 


were rounded off to 5 decimal places prior to output. 


3. Least Squares Adjustment by the Condition Equation 


— = — —  — шшш» ee — ee ee ee — ww SE a no шш 


The general principle of least squares adjustment by 
the condition equation method in surveying is 
onem еа оошсо Сота уос о the sam of the 
Squares of the corrections to the observations plus tke 
necessary mathematical conditions involving some or all 
MremGOrrect1Ons; each condition by itself is made ecual 
Me ro bY adg@ing the Corrections ts the discrepancy 
determined from a preadjusted calculation (i.e. caiculated 
ШОО e Observed values, rather tnan from adjusted | 
values); thus, the magnitude of the sum of the squares is 
meoewechangead by adding Conditions [Ref. 8]. 
For uncorrelated observations of unequal precision, 
mie principle condition SZ least squares function пау Бе 


expressed in matrix form as follows: 
0 = vVipvy - 2kK(BV * RF] 
= vVipy - 2[BV + W]'K 
= vipv - 2¢(Bvy’ + wl yK 
= vipy - 2[VTBT+ gw! Jk 


=i ep awe ci 


minimum 


where U is the least squares function matrix, K is a matrix 
Of Lagrange multipliers, P is the weight matrix, V is the 
wee on vector, B is the constant coefficients of the 
Corrections, and W is the precalculated discrepancy 

NIE. 9]. It is numerically more convenient for later 
Seyelopment to multiply by 2. Taking the partial 
derivatives of the U matrix with respect to each of the 


corrections and equating to zero leads to 


4 7 


сә 
c 
II 


22 Ve een = 0 


Q 
< 


OE 


ÇV = Р-1В К (3.26) 


This equation is called a correlates ecCtlati on orecor ea 
equation matrix. The solution of the Lagrange ов ар 
Matrix can be derived by multiplying Equation (3:26) Бу 


B matrix then adding the W matrix 
BV + # = ВР-1В К + к 


(BV + W) is the condition equation matrix which must equali 


zero. The solution of the Lagrange nultiplier vectores 


BP -1BÍ K 


=N 


K (BP-ipl) -1 (-wW) (3. 27) 


The correction vector is derived by substitutimsa 


Equation (3.27) into Fauataon nee) 
V = P-1Bl(BP-187)-1 (- x) (3.28) 


For uncorrelated observations of equal precision, 
the correction vector can be derived similarly to the 
unequal precision case or by uSing the inverse of the weight 
matrix which is the identity Matrix. qun ion m p 


becomes 
V - B! (BB ) -1 (- x) (3.29) 


Closed traverse station positions may be adjusted bv 
using the technique of least squares adjustment by the 
Condition Equation Method. There are two condition equa- 
tions. They are the azimuth and coordinate condition equa. 
tions. The coordinate condition ejuations are divided ies 


two parts, which are E and N coordinate condition есоа стони 
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The azimuth condition equation is the sum of the 
eorr ections O the angles in a traverse plus a precalculated 
discrepancy W, (Equation 3.6) which must equal zero. Where 
О К едпа1< the corrections to the angles, the Jeneral equa- 


at 
tion can be written: 


РЕ n n-1 
ОУ eee + W. = O (3-30) 
i - 1 1-1 

where V3; is the correction to the distances. However, 


this term equals zero and does not effect the equation, but 
it does reserve space in matrix notation for distance 
corrections. A closed traverse which was conducted at “oss 


ending can be written with an aziauth condition of 


à + à + = + ы + с + z = 
1 ү 1 У 1 ДЕ 1 m 0 E 0 Yao + 0 V44 


+ W, = 0 (3-31) 


Meee the precalculated discrepancy (f4) is +8.844" or 
+0.00004288 radians. 

The Simple linearized form of an E coordinate condi- 
tion equation is the sum of the adjusted departures and must 
БО zero. sy applying Equations (3.2) and (3.7), the 


emal formula for E coordinate condition equation is 


— 1 


where : 
ON Rc EZ + Lia, Еи вов (3.33) 


ai 
= 1] 

fees the adjusted azimuth. By using Equations (3.32) and 

(ees) With the data conducted at Moss Landing, these 

equations become 


(di, + v,,)-sin Az., + (4, + V42)-Sin Az,, 


* (d. + V 3). Sin Az. БИЕ 


А уз (3.34) 


T 


49 


where 


c ~ ВЕЕ š A, + Vat 
= AZ, + Va, 
AZo = Age + © е оону MIU 
- A2 t TENE 
AZo, = Aze + A, + v, * C, * v5 * € .4* v s 
= aZ. t REL NE. 
Az. is the precalculated azimuth or unadjusted 
az1müth (Eguation 596 


From trigonometry 


sin ( А + (А ) 
іп which /NA is a very small angle in radians, 


cos/NA = 1 land) sumo A 


therefore, 


sin ( A */NA ) 


sino ce d AC SEMEN 


fcon the first tern o: Equation ao) 
СО 

E as vj.)sin (Azo ey) 

5 у 41) 5220 ÀZ, + у Соз Àz,] 

= d,sin Az, + v,, d,cos Az, + v,,Sin Az, + VY, 


Ср is very small, by letting this Pern 


Equation (3.35) becomes 


(Чү * v4,)Slin Az, 


] + + 
d, sin Az, v. , d, cos А2, ү 


41 91h AZ, 
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sin À.cos/NàÀ + cos A.sin/f\A 


let 


cos AZ, 
(3.355 


equal Zener 


Уз 1 


(3.36) 


пле сесола апа ЕҺі Еа terns of Equation (3.34) сап be derived 


Еш ат tO Equation (3.36), so 


(15 + узо) іп А2. 


= d, sin А2, + (V. * y,5,)05cos А2. t v o Sin AZ, ОЗ 
(d. * V4) Sin Az., 
= d, sin Az, + (V. i *V,5 *V, 4) d4cos AZ, + Vj, Sin AZ, (2. 39) 


“i SP nh in mF ati onmsm( 3755), (3.37), and (3.3%) 


mimeo Equation (3.34) and rearranging terns 
V41 (d1cos Az, * d5cos Az, * d4cos Az.) 
t v» (d „CoS Az, + d4cos А2.) + vV,404Ccos Az, 


* V Sin AZ, + v Sin AZ, + v 


di 42 Sin А2. 


аз 


* d,sin Az, + d,sin Az, + d,sin Az, - (E;- Еу) = 0 


This equation can be written in the gereral formula as 


=! n-1 n -1 n -1 
Eu. Pos Az.) + 23 San А2. + > asin AZ. 
izl К= і p e E | 

EE E = 0 (3.39) 
ER ting socunatlouso(3-2), (3.3), and (3.7) into Equatior 
(2-39) 


n-1 


p! 2! 
2.08 251 E АШ и = 0 (3.40) 
= = = 


For example, if the number of the observed angles (гп) at 


Moss Landing is four, this equation can be expanded to 
+ 
(AN, +N + SR.) NEN s N HON) V o 


* 1 + l 1 = 
v44 Sin AZ, Vio Sin AZ, + Vj3 SLD Az. * Wo 0 
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Substituting the numerical values from the precaiculated of 


this traverse into this equation, the result is 


(6702.37612)v,, +(5318.48345)v,5 +(4999.28284)v¥, + O.¥, 


-(0.23567) у, *(0.88818)v., *(0.63570)v,, * 8, —- 0 (3.41) 


where the precalculo t oS SC repe (84) is 0.08608 п аве 
0.v.,is equal to zero, which does not effect this equation, 
but space is reserved in matrix notation for the last angle 


cCcorrection, 


The N coordinate condition equation is the sum o£ 
the adjusted latitudes and must equal zero. By applying 


Equations (3.3) and (3.7), the general condition equation 


- N.) 00 


-1 
(d. + ху) COS А2 4: - (NIE 
T 


This equation can be derived similarly to the E cooramwmr 


condition equation сасе сезше иаи 
nex 


noy noe : 
-2( DAE) * Dv cos Az; *W, -0 (3.42) 
е] Kod = 

From the data acquired at Moss Landing, this equation can 


Le expanded to 
та (E, + OE, + OE) 7 Vào (ES nom E OF, 


+ 2 + + + Wa = 
za COS À 1 E cos AZ» з 69S AZ. 3 0 


Substituting the numerical values from the calculation еб 
this traverse into this equation, the result is 

~ (3959.89460) v,, ~(4295.49626) v_. 7 (9116294733) ME + 0- v, 
*(0.97183) v. + (0.87274) у, * (0.77194) v. + K. = (3. 43) 


where the precalculated discrepancy (8W4) is -0.089236 m. 


3 


ЕКОЛО О апа (3.42) can be written 
im the matrix form as (BV + W), where B is the constant 
БЕБИ Са еее ое Еле corrections, V is the correction vector 
which is composed of the angle and distance corrections, FH 
is the precalculated discrepancy. Equations (3.31), (3.41), 


and (3.43) are written in matrix form as 


det 02207612 3959559460 V. 
1 5318.48345 -4295.49626 Ves 
1 4999.28284 -4116.94755 V. 
BT = 1 0.00000 0.00000 v = V. 2 
0 -0.23567 0.97183 v 
0 0.48818 0.87274 Vy 
0 0.63570 0.77194 v 
-W, = -0.00004 
Шы = | -я„ = -0.оввоа 
-# = 0.08936 


For uncorrelated observations of unequal precision, 
the inverse of the weight matrix is the diagonal matrix, 
which is equal to the observation's variance. The diagonal 


elenents of the inverse of the weight matrix is 


р-1 = sin2(1.984") = 0.92519 х 10-19 
p^, = зіп? (1.405") = 0.46398 х 10-109 
pot = sin?(1.203") - 0.34016 x 10-19 
р-1 = sin2(1.614") = 0.61229 x 10-10 
ОК = 0.10000 х 10-5 
Е ООШ 92 00000 0105 
RC - 0.90000 x 10-5 
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For uncorrelated observations of unequal precision, 


the correction vector is found by solving EQuaGrOn ЕУ 


Vi, = ee eS a ee radians or 12920090 

ym -1.631140 x 10-5 radians or -3.36446" 

Уз “тозо NE БЕ radians or =2. 647908 
ү = V. 4 c2: radians ос 4.7531 

Oro: 0.00924 meters 

“чо = 0.00039 meters 

V44 0.00357 meters 


The corrected angles are 


C = 2460 05' 45. 12117" 
An = 222095 1 05 o S 
CA, = 1909 14' 59.95242" 
CL, = 2779 05' 12.24688" 


The Corrected distances are 


d. = 1424.00424 ao. 
dj = 365.74439 m. 
da = 6476.27457 а. 


Using the corrected data to recompute the cooriinates is 
similar to the initial traverse computations. The adjustegs 


cCosordinates are 


Eo = 607,943.45521 m. 
н, = 4,073,939.74809 m. 
E, = бос оте 
N, = 4,074,258.95029 m. 
E, = 612,238.85256 m. 
н, = 5,079,258.31758 ш. 


The last adjusted coordinates are the same values as the 
fixed coordinates. Hence, it is not necessary to iterate. 
For uncorrelated observations with egual preci- von ee 


correction vector is found by solving Equations: 
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The corrected angl 


90292 


М u N 


e 


1.26884 
230226550 
- 1. 751808 
- 223 00 

0.14034 

0.23842 

0520362 


are 


2460 0 
222975 


The corrected distances are 


d 


1 


5 


d 


By using the corrected 


É 


the adjusted coordinates are 


52 
Ho 
E3 
Ns 
E4 
NA 


b P? b b bd b b 
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1° 


1900 14' 
2779 05! 


КОБЕ ааа 
10-5 radians 
10-5 radians 
10-5 radians 
10-7 meters 
10-7 meters 
107" meters 
45. 81717" 
05.36599" 
5991906" 
1225 9250." 


T424.00400 n. 
365.74400 n. 
6476.27100 m. 


data to recompute 


607,943.45994 m. 
20g 93914899 m. 
50:8 012 220 0T mn. 
47074, 256295013" m. 
ОИС 2 6165256” п. 


= ü,079,258.3175% n. 


OL 


Or 


OF 


Oc 


geo ile 17" 
512 2501 
5-601909 
E5075 


the cooriinates, 


The last adjusted coordinates are the same values as the 


med coordinates. 


Hence, it is not necessary to iterate. 
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The standard deviation of an observation which has 
unit weight can be fend, спе гоо equations. Mom 


uncorrelated observations of unequal precision, 


6, - wey) 173 (3.44) 


r 


and, for uncorrelated observations of equal precision, 


б. = Бар (3.45) 
r 


where r is the number of condition equations. There are 
three condition equations for a closed traverse (r53) g 
standard deviations of adjusted angles and distances can be 


found by the following equation: 


6. = T 172 (3.46) 


where 6. is the standard deviation of the ith adjusted 
quantity. The quantity in the ith row of the оос е Qu 
is an elenent of the P-! - p-1gl (ge-1g!) -1gp-1 matrix rom 
uncorrelated observations of unequal precision 

(Ref. 3, p. 918]. For uncorrelated observations Of yeaa 
Precision mon is the element of the I - B! (BBĪ) -1B matrixa 
The matrix Р-1 - р-1В!(Вр-1ВТ)-18р-1 ог ёһе matrix 

I- B'(BBT)-1B is equal to 


Q 44 012 oy See ON. Qin 
“a U a e 
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For uncorrelated observations of unequal precision, 


mae Standard deviation of unit weight is round by solving 


Equation (3.44), the result is 


SEC 


° = 


2269679 


The results of P-! - p-18! (Bp-13T)-15p-1 are 


ОЧИ ОЕ Е ОООО ОО 5. T595-—6-963 -63.583 

ЖО Оо о UMS 0.902 32915 85.6989 41.59 

B ovd» 0.29]05-0.0149 2.5969. 3.9578 31.79 

10—109 ОСЕТИТИ 00 0 0 11950.0027 1.513 1.3501 -9.500 
ОИ m» 5563951.913—— 99172 —117.8 -1097 
oO. Oem S Eoo 0T TTL. 9830 -1525 

SG Sep с д | ons 0 0500 — 15403. —1525 152 3 


The standard deviations of adjusted angles and distances are 


pound by solving Equation (3.56); the results are 


Os, = £ 1.803400 seconds 
68 = + 2.64096 seconds 
O33 ani 2.59253% Seconds 
om = + 0.26136 seconds 
Og, = + 0.00269 meters 
СЕА ООСО ОС ЕЕЗ 
баз = + 0.00745 meters 


For uncorrelated observations of equal precision, 


the 


Seamdard deviation of unit weight is found by solving 


Eaton (3.585), 


б 


о 


the result is 


2 


а Е в 


The results of I - B'(BBlj—1!B are 


1.85972 -27.199 26.962 -1.235 -0.5T70N DM M UM 
-27.199 503.1 -498.7 22.855 ОЗО DL SE 
26.962 -498.7 494.4 -22.66 01570 O K 0и 


10-3 -1.235 22.855 -22.66 1.036 -9: 20 P s eee 
zs 0.348 0.370 -0.201 10000-0001 EDD 
07903 0.880 0.523 -0.159 -0.001 ОЕ E 
205955 0.473 0.518 -0.134 -0.001 02001 itis 


The standard deviations of adjusted angles and distances are 


found by solving Equation (3.46); the results are 
бат = + 0.15918 seconds 
On = + 2.94262 seconds 
e = + 2.91722 seconds 
Ӧд = + 0.13369 seconds 
Qr. = + 0.00002 meters 
Oyo = + 0.00002 meters 
баз = + 0.00002 meters 


The computations in this subsection were also 


performed on the NPS IBM 3033 computer; 


places, 


finat solutions. 


they used 165 dcecimal 
and were rounded off to 5 decimal places in the 


These are more decimal places than are 


normally used in practice because such precisions are not 


generally attainable Ey corresponding observations. 


However, in comparing two computational methods tatg 


theoretically equal, 


rounc ott error than the other 


in the next chapter. 


5B 


one method may be more sensitive to 


This will be commented on 


IV. DISCUSSIONS AND ANALYSIS OF RESULTS 

Three programs were used to compute the initial traverse 
em@dead just the closed traverse station positions. Ргодгапѕ 
1, 2, and з were used to compute and adjust the traverse 
eae ton positions by the Approximate Method, the Indirect 
Observations Method, and the Condition Equation Method, 
respectively. The programs were written in WATFIV language 
meena telementation on the NPS IBM 3033 computer. The 
Kas er output and listings of Programs 1, 2, and 3 are 
provided in Appendices A, B, and C, respectively. The 
maximum number of intermediate traverse station positions 
that can be computed and adjusted by these programs is 30. 
The programs were tested using several fictitious data sets 
to ensure their performance in handling the various 
intermediate traverse station positions. 

The computer storage area and CPU time of Programs 1, 2, 
and 3 has been compared (Tabie X). Тһе adjustment of a 
Closed traverse by the Approximate Method did not use a 
weight matrix in the matrix computations.  Tne computer 
program for this method was written by using the variables 
in only one dimension for storage of both data and results. 
This method did not require any iterations. The least 
Squares adjustment of a closed traverse characteristically 
is used to simultaneously eliminate closing errors in 
azimuths and distances, and Programs 2 and 3 both utilized 
such an adjustment. The computation by either Least Squares 
Method used matrices in two dimensions to compute the 
correction vector. Likewise the computer programs were 
written using variables in two dimensions. Included are all 
subroutines from Program 1 and extra subroutines for each 
individual program. Therefore, both Programs 2 and 3 used 


ШШ FS Computer storage area and CPU time than Program 1. 
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The Indirect Observations Method requires that the 
number of observation equations be equal to the sum ot 
observation equations of both observed angles and distances. 
The number of observation equations is not constant. It 
changes depending on the number of the traverse stations-- 
Bags Means that the row dimension of each matrix wili change 
depending on the number of observation equations. 

Therefore, the subroutines affected by the number of trav- 
erse stations are difficult to write in the general program. 
The coefficients in the A matrix are computed by taking. the 
partial derivatives of the functions with respect to each 
unknown variabie, where the number of unknowr variables is 
not constant. The number of unknown variables changes 
depending on the nunber of the traverse stations. 

Therefore, the subroutine for the computed A matrix is also 
ЕЕЕ TC CU]t to write in the general program. The Condition 
Eguations Method has only three equations which makes it 
easier to derive the general forn and write the computer 
program. The adjustment using the Indirect Observations 
Method does not apply the corrections to the observed values 
directly; therefore, it requires the approximate values of 
the unknown coordinates be computed for the correction 
vector of the unknown coordinates. This method requires at 
least two iterations to check the insignificance of the 

enr ection vector when it is compared with an arbitrarily 
selected small number. The adjusted coordinates from the 
Condition Equations Method can be checked at the first iter- 
BENOIT Thus, Program 3 is more economical as it uses less 
computer storage area and CPU tine than Program 2 (Table X). 

The closed traverse at Moss Landinj originated and 
terminated at control points with known positions which were 
determined by third-order methods. The azimuth closure was 
EEUU" per station and the position closure was 1:66,617. 


E eclassritication, standards of accuracy, ani general 
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specifications for a third-order class I traverse s m 
indicate the azimuth closure is not to exceed 3" per station 
and position closure must be Letter than 1:702 000 s 
traverse met the specifications and standards of accuracy 
for a third-order class I traverse. 

Using Programs 1, 2, and 3, the traverse was computed 
and adjusted in UTM grid coordinates. The finals G 
Coordinates were transformed from UTM grid coordinates to 
geographic coordinates [Ref. 4, рр. 319-321]. шеш еы 
ence between the coordinates for each nethod and NOS 
[Ref. 10] have been computed (Table XI). The technique of 
Least Squares for both methods yields identical computa- 
tional resuits to at least five decimal places (Table Xi, 
Methods 3.1 and 4.1; 3.2 and 4.2). Least Squares proves 
the Lest estimates for positions of all traverse stations. 
Program 2 performs a statistical test, yielding the standard 
deviations of both adjusted positions and observed values. 
Conversely, Program 3 only yields the standard deviation of 
adjusted observed vaiues.- Comparisons of the standard devi- 
ations of adjusted observed angles and distances were made 
between these computed by the Indirect Observations meas 
and the Concition Equations Method (Table XII). The 
computed standard deviations of the adjusted angles ditfered 
significantly in the fourth decimal place in several cass 
which is a slight indication that the Indirect ObServabusm 
Method is more sensitive to round off error. The standard 
deviation estimates are larger for the Indirect Observations 
Method in every case of a significant difference in the 
fourth decimal place. However, a significant difference in 
the fourth decimal place is insignificant in terms of the 
okservational precision. The standard deviations of the 
adjusted distances for both methods yield identical resuits 
(Table XII). Both Methods differ from the NOS positions 


the fourth decimal place in seconds of arc for both latitude 
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EN lon eaen C UTheShorrzontal control for third-order stan- 
dards recuires accuracy to three decimal places in seconds 
Beyare for both latitude and longitude--thus, these methods 


can be used for computing third-order positions. 
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АС NSGONCTUSEONS 


The conputer programs developed by the author ard 
contained ir this.thesis have a variety of usəft31 ani eE 
tical hvéárographic applications. In hydrography geom 
rilelG work plays a key role.  dorizontal control is ness 
Sary in determining positions andi, even in areas which 
appear to initially have adecuate control; au UON M 


control stations often need to te established alitcr o ee 


in tre 2ield. These field-genrerated control points SUUM 


in 


adjusted to Zit within the existing] sarye? net an: SE 
tize errors iu tne 2lell Aeasurescrts. “The errors cana 
minimized by these computer programs. Direct application of 
these programs would be enormcusly beneficial, since no 
equivaient software exist at MPS to perform such adjustments 
at the present time. 

In the Approximate Method, a traverse is adjusted by 
computations using a hand calculator. This method iS MEM 
able for field computation and the adjusted coordinates meet 
the specifications and standards of accuracy гог а бои 
order class I traverse. Although the accuracy rong ss 
method is less than the Least Squares Metnod, it does not 
require a computer and it furnishes coordinates which@e@auee 
checked in the field by the fred Darty: 

The Least Squares Method provides the best estimates for 
rositions of all traverse stations. Least Squares does 
require pore of a computational effort than the Appro mick 
Method. However, the accuracy of Least Squares is better 
than the Approximate Method, and the Least Sguares Method 


should be used for the final office Computations ME 
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Least Squares Methods yield identical results, but the 
Condition Equations Method is more economical and easier to 
derive than the Indirect Observations “Method. Therefore, 


meme COndition EGUations Hethod should be applied at NPS. 


B. RECONMENDATION 


Studies should be continued at NPS to compare the 
economy of adjustment methods used in this thesis with other 


methods or least squares adjustment. 
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